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RESUMEN 

Se han seleccionado cinco sismos interplaca tipo thrust, los mismos que 
han sido escalados al espectro de la Norma Ecuatoriana de la Construccion de 
2015, que corresponde al sismo de diseno DBE y al espectro maximo considerado 
MCE; para encontrar la respuesta en el tiempo, mediante analisis no lineal de cada 
uno de los porticos del Bloque Estructural 3. 

Se seleccionaron cuatro sismos ocurridos en Peru y son los registrados en 
1966, 1970, 1974 y 2007. El quinto corresponde al sismo de Bahia de Caraquez 
de 1998, en estos sismos se aplico la tecnica “weighted average” para el 
escalamiento de acuerdo a lo que propone el ASCE 7-10. 
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Estos sismos son de larga duracion, con una o dos fases intensas y con 
varios pulsos fuertes de energia que los convierten en sismos muy destructives, 
por lo que era necesario conocer las respuestas en aceleraciones y 
desplazamientos del aislador y de cada uno de los pisos de la estructura. Por otra 
parte se determino espectros de piso que serviran para estudios futuros de calculo 
de perdidas. 

Palabras claves: Escalamiento de sismos. Respuesta en el tiempo. Metodo 
“weighted average”. 


ABSTRACT 

Five interplate thrust earthquakes were selected, The seism were scaled to 
the spectrum of the "Norma Ecuatoriana de la Construccion 2015”, that 
corresponds to the design seism, DBE and to the maximum considered seism, 
MCE; to find the time response using inelastic analysis of each frames structural 
block 3. 


We selected four seism of Peru that were registered in 1966, 1970, 1974 
and 2007. The fifth seism corresponds to Bahia de Caraquez of 1998, for the 
scaling of these seism we apply the “weighted average” method, recommended by 
ASCE 7-10. 

These earthquakes are long period, with one or two intense times and 
some strong pulses of energy that convert in seism more destructive, for this 
reason we need to know the accelerations and displacements of the insulator and 
each floor of the structure. On the other hand we found Spectrum of floor that will 
be using to future studies of calculating losses. 

Keywords: Scaling of earthquakes. Time response. Weighted average method. 


1 INTRODUCCION 

En terminos sencillos se puede decir que los aisladores sismicos 
desacoplan a la estructura del suelo ya que se crea una interfaz entre estas. 
Interfaz creada por los aisladores que tienen una rigidez lateral muy pequena que 
ocasiona una flexibilizacion de la estructura y por ende un mayor periodo natural 
de vibracion. 

Al tener un mayor periodo de vibracion se ingresa a una zona del espectro 
donde las ordenas espectrales son muy bajas pero las ordenadas del espectro de 
desplazamiento son muy altas. 

Para ilustrar lo expuesto en la figura 1 se presentan los espectros de 
aceleraciones, a la izquierda y de desplazamientos a la derecha para un Proyecto 
Estructural ubicado en la ciudad de Quito, el espectro de aceleraciones se obtuvo 
utilizando la Norma Ecuatoriana de la Construccion, NEC-15, en un perfil “D”. En 
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cambio el espectro de desplazamientos se obtuvo de un estudio de peligrosidad 
sfsmica para ese sitio, en que se modified la forma del espectro de 
desplazamientos. 
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Figura 1 Espectro de Aceleraciones y de Desplazamientos, para un Proyecto 

Estructural de la Ciudad de Quito. Fuente: Aguiar y Morales (2015) 

Lo que se desea ilustrar es que si el periodo de vibracion de la estructura 
con aisladores es 3 segundos, por ejemplo, la aceleracion espectral es baja, 
alrededor de 2.5 m/s 2 pero el desplazamiento espectral es aproximadamente de 
48 cm. De tal manera que las fuerzas son bajas pero los desplazamientos laterales 
son altos, en el sistema de aislacion, lo que debe ser cubierto con un aislador que 
tenga una capacidad de desplazamiento mayor a la demanda. 

Los aisladores sfsmicos incorporan amortiguamiento a la estructura, en 
una cantidad considerable, claro esta dependiendo de las caracterfsticas y tipo de 
aislador, tema que se hablara posteriormente. Por ahora, lo que se desea dejar 
claro es que este amortiguamiento de los aisladores conduce a menores 
ordenadas tanto en el espectro de aceleraciones como en el espectro de 
velocidades. Aguiar etal. (2008). 




Entonces, lo que preocupa en las estructuras con aisladores es conocer 
en primer lugar la respuesta que van a tener ante la accion de diferentes sismos y 
segundo si el aislador seleccionado tiene la capacidad de desplazarse lo que la 
demanda sfsmica le impone. 

De ahf la importancia de realizar estudios de peligrosidad sfsmica de un 
sitio del proyecto, para definir sismos con fuente similar, y encontrar la respuesta 
no lineal de la estructura a sabiendas que son los aisladores los que van a 
incursionar en el rango no lineal, tema que sera abordado en los siguientes 
apartados de este artfculo. 

Es importante destacar que a nivel mundial, existe conciencia de que la 
mejor forma de disenar las estructuras es mediante la incorporacion de sistemas 
de control pasivo a la estructura. Sean estos aisladores sfsmicos como en el 
presente caso (Retamales et al. 2015; Constantinou et al. 2011; Sarkisian et al. 
2011). O mediante la incorporacion de disipadores de energfa. (Bozzo y Barbat, 
1999; Hurtado y Bozzo, 2008; Nuzzo etal. 2015). 
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2 SISMOS INTERPLACA TIPO THRUST 

El ultimo sismo fuerte, interplaca tipo thrust (en ingles) o con bajo angulo 
de buzamiento, con magnitud mayor a 7, registrado en el Ecuador es el de Bahia 
de Caraquez del 4 de agosto de 1998, que tuvo una magnitud de 7.2 y una 
profundidad focal de 22 km. La Intensidad en la Escala de Mercalli Modificada de 
este sismo en la zona epicentral fue de VIII, como se observa en la figura 2; la 
isosista de intensidad VII cubrio una buena parte de Manabf; la curva de 
intensidad VI llego hasta Guayaquil, que se encuentra aproximadamente a 250 
km., del epicentro. 



Figura 2 Isosistas del sismo del 4 de agosto de 1998 y replicas hasta diciembre de 

1998. Fuente: Egred (2009). 

Lo que se quiere destacar con lo expuesto en el parrafo anterior, es que 
los sismos interplaca tipo thrust no solo causan dano en la region epicentral, sino 
que tambien producen dano a grandes distancias. 

El sismo de Guayaquil de 1942, es otro sismo interplaca tipo thrust, su 
epicentro fue en Manabf y tuvo una magnitud de 7.9. Se conoce con el nombre de 
una ciudad distante del epicentro por el gran dano que causo en las pocas 
construcciones de hormigon armado de esa epoca, a tal punto que se llego a 
pensar que este tipo de material no es apto para zonas sfsmicas. 

El sismo de Maule del 27 de febrero de 2010, que tuvo una magnitud de 
8.8, es otro ejemplo de que el dano no solo se produce en la zona epicentral sino 
que alcanza grandes distancias. En efecto este sismo llevo al colapso el edificio de 
altura Alto Rfo ubicado en Concepcion y genero mucho dano en otros edificios 
altos de Santiago, ubicados a 400 km., del epicentro. (Aguiar 2010; FCFM 2012). 

Los sismos interplaca tipo thrust de magnitud mayor a 7, se caracterizan 

por: 


Ser sismos de gran duration. 

Generan dano a grandes distancias del epicentro. 
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• Tienen varias fases de movimiento intenso. 

• Sus espectros en la zona descendente tienen varios picos caracterfsticos, 

en los que se incrementa la aceleracion espectral. 

En el Ecuador se tienen cuatro registros del Sismo de Bahia de Caraquez 
de 1998, dos de ellos fueron registrados en Manta y Salinas pero tuvieron 
problemas y no sirven ya que llegan a un valor maximo de 0.6 g\ no es que este 
sismo tuvo aceleraciones espectrales mayores al 60% de la gravedad, sino que 
hubo algun desperfecto en el registro. 

Otro de los registros fue obtenido por la red sfsmica del Austro en Cuenca, 
pero tiene valores muy bajos de aceleracion. Finalmente el ultimo registro y es el 
que se utiliza en este estudio fue hallado en la Universidad Catolica Santiago de 
Guayaquil por el Instituto de Investigaciones de la Facultad de Ingenierfa Civil. Se 
tienen las dos componentes horizontales que tuvieron aceleraciones maximas de 
120.9 gals, en la componente Norte Sur y 191.9 gals en a componente Este 
Oeste. La componente vertical tuvo una aceleracion maxima de 85.6 gals. La 
duracion de estos registros fue de 90 segundos. 

En este artfculo, interesa encontrar la respuesta en el tiempo del Bloque 
Estructural 3 de la Universidad de Fuerzas Armadas ESPE, ante algunos sismos 
interplaca tipo thrust cuya fuente tectonica sea similar a la de Ecuador y lo mas 
apropiado para este efecto es trabajar con los sismos registrados en el Peru y que 
se indican en la figura 3. Por dos motivos, el primero de ellos por su cercanfa con 
el Ecuador y el segundo porque fueron de magnitudes entre 7.9 y 8.1, similares a 
las que se esperan en el Ecuador para el sismo de diseno. DBE. 
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Figura 3 Sismos Interplaca Tipo Thrust seleccionados en este estudio. 
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En la tabla 1 se presenta informacion sismologica basics de los sismos 
seleccionados; las zonas afectadas y el numero de vfctimas de estos terremotos. 
Los sismos de Peru de 1966, 1970 y 1974 fueron utilizados en la elaboracion de la 
norma sfsmica de Peru de 1997. Completa la lista de sismos seleccionados el de 
Pisco de 2007, que acabo con muchas construcciones de adobe. 


Tabla 1 Ano, ubicacion, magnitud, profundidad focal y afectacion de los sismos 
de 1998, 1966, 1970, 1974 y 2007 


ANO 

FECHA 

COORDENADAS 
DEL EPICENTRO 

MAGNITUD 

K 

PROFUNDIDAD 

FOCAL 

ZONAS AFECTADAS 

VICTIMAS 

ECUADOR- 

1998 

4 de 
agosto 

0° 33' 0" S, 
80° 31' 0" W 

7.1 

22 

Bahia de Caraquez, 
Canoa; San Vicente 

1 muerto, 100 heridos 

PERU-1966 

17 de 
octubre 

10° T 0" S, 
7° O' 0" W 

8.1 

37 

Costa del 
departamento de 
Lima, entre Lima y 
Supe 

220 muertos, 1.800 
heridos 

PERU- 1970 

31 de 
mayo 

go 2g , 12 „ ^ 

78° 18' 36" W 

7.9 

64 

Ancash, La Libertad, 
Huanuco y Lima 

70 000 muertos, 50 
000 muertos y 20 000 
desaparecidos 

PERU-1974 

3 de 
octubre 

12° 15' 0" S, 
77° 31' 12" W 

8 

21 

LIMA, CALLAO, 
COSTA CENTRAL 

252 muertos, 3600 
heridos 

PERU-2007 

15 de 
agosto 

13° 21' 10.8" S, 
76° 30' 36" W 

7.9 

40 

Costa central del 
Peru 

595 muertos, 2 291 
heridos, 431 000 
damnificados 


A continuation en la figura 4 se presentan los acelerogramas de los 
sismos seleccionados, a la izquierda se tiene la componente Este-Oeste y a la 
derecha la componente Norte-Sur. Se aprecia que son sismos de larga duracion. 
En el de Bahia de Caraquez de 1998, en el de Callao de 1974 y en el de Pisco de 
2007 se observa en el registro, como que se habrian dado dos sismos en uno 
solo; en efecto se ve que una vez que termina la fase intensa y parece que el 
sismo finaliza, nuevamente se activa con otra fase intensa. 

En todos los sismos seleccionados la fase intensa es de larga duracion. El 
sismo de Pisco de 2007 duro alrededor de 3 minutos. Por cierto, su registro fue 
obtenido en la Universidad San Luis Gonzaga de la ciudad de lea, situada 
aproximadamente a 240 km, del epicentro. 

Interesa conocer cual va a ser la respuesta de la estructura del Bloque 3 
ante estos sismos que son de larga duracion, con fases intensas muy grandes y 
que son dos sismos en uno. Interesa ver la respuesta no lineal. 

En la tabla 2 se indica la nomenclatura con la cual se va a identificar cada 
uno de los registros, las aceleraciones maximas y minimas y el tiempo de 
duracion. Las aceleraciones son bajas porque los registros fueron obtenidos a 
distancias grandes de la fuente y para utilizar en la evaluation de la respuesta no 
lineal del Bloque Estructural 3, de la Universidad de Fuerzas Armadas ESPE, 
tienen que ser escalonadas (incrementadas) para que sean mayores que las 
ordenadas espectrales de la Norma Ecuatoriana de la Construction, NEC-15, que 
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es el sismo DBE y mayor que las ordenadas espectrales del sismo MCE que se 
obtuvo multiplicando el espectro del DBE por 1 .4 


BUM.n MKA.IO 



t St M M H 0 It St to M M 

TM VM 




Figura 4 Acelerogramas de sismos interplaca tipo thrust seleccionados para 
encontrar la respuesta no lineal de las Estructuras del Bloque 3. 



Roberto Aguiar, Majorie Arias, Gregory Pachacama y Evelyn Iza 


158 


Tabla 2 Aceleracion maximas y mfnimas de los registros considerados en el 
estudio y nomenclatura que se utiliza en el artfculo. 


Registro 

Ano 

A 

** max 

(gals) 

A . 

** min 

(gals) 

Tiempo 

(seg.) 

Bahfa_NS-88 

1998 

121 

115 


Bahfa_EO-88 

1998 

192 

115 


7035-66 

1966 

168 

181 

65,66 

7036-66 

1966 

247 

269 

65,66 

7038-70 

1970 

93,8 

105 

45,12 

7039-70 

1970 

91,2 

97,8 

45,12 

7050-74 

1974 

179 

127 

97,86 

7051-74 

1974 

159 

192 

97,86 

7081-07 

2007 

272 

260 

163,68 

7082-07 

2007 

334 

301 

163,68 


3 ESCALAMIENTO DE SISMOS 

En la figura 5 se presentan los espectros de respuesta elastica para 5% de 
amortiguamiento de los sismos seleccionados; para cada una de sus componentes 
horizontales de movimiento del suelo. Se aprecia que todos ellos tienen varios 
picos caracterfsticos en la zona descendente del espectro y el de Peru de 2007 
tiene valores altos en la zona de periodos comprendidos entre 2 y 4 segundos; 
este comportamiento no consideran la mayor parte de normativas sfsmicas por 
esto es que Saragoni y Ruiz (2012) recomienda tener espectros con dos 
plataformas de aceleracion maxima para los sismos interplaca tipo thrust. 

En la parte inferior de la figura 5, se presentan todos los espectros de 
respuesta y se ha encontrado el espectro promedio, el mismo que se compara en 
la figura 6 con el espectro que se halla al aplicar el NEC-15 en un suelo tipo C 
(Espectro de diseno). 

Se aprecia que el espectro promedio de los 10 registros considerados se 
halla bastante distante (menor) que el espectro del NEC-15 por lo que los registros 
tienen que ser escalados para que sean mayor o igual al espectro del NEC-15 
pero no en todo el rango de periodos sino en el rango que recomienda el ASCE 7- 
10 . 

En efecto, el ASCE 7-10 recomienda que para estructuras con aisladores 
sfsmicos el rango de interes se encuentra entre 0,5 T DBE < T < 1,25 T MCE \ siendo 
Tdbe el periodo del sismo de diseno y Tmce el periodo del sismo maxi mo 
considerado. 

El periodo de vibracion del Bloque Estructural 3 que tiene aisladores 
sfsmicos es T = 3.5393 s para el sismo DEB y para el sismo MCE es T = 3.6263 s. 
Por lo tanto el rango de interes para el escalamiento es: 1.769 s <T < 4.532 s . 
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Figura 5 Espectros de respuesta elastica de los sismos considerados y 

Espectro promedio. 
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1.40 



T(s) 


Figura 6 Promedio de los espectros de respuesta elastica de los sismos 
considerados y espectro de la Norma Ecuatoriana NEC-15, expresados como una 
fraccion de la gravedad (espectro de diseno DBE). 


Para minimizar la diferencia que existe entre espectro objetivo A (espectro 
de diseno del NEC-15 para el sismo DBE) y el espectro de respuesta A se calcula 
el factor Fj. La manera como se obtiene es la que da origen a los metodos de 
escalamiento. Asf el ASCE/SEI 7, propone que se obtenga con la siguiente 
ecuacion. Reyes y Kalkan (2012) 

Ej = min F1 ||log(i) - log(F ; - ^)|| (1 ) 

Donde || . || es la norma Euclidiana. Se requiere para este proposito un 
metodo numerico que minimice las funciones escalares de una variable. La 
minimizacion asegura que cada espectro de respuesta escalado se encuentra muy 
cerca del espectro objetivo, como se muestra a la izquierda de la figura 7. 

Un procedimiento mas especffico para calcular los factores F,, que se 
utilizara en este artfculo es el denominado “weighted average” en este caso se 
tiene: 


Ej =J> [Mt,) - Fj AiT,)} 2 (2) 

i = 1 

Donde son pesos que se asignan para los perfodos seleccionados 7); 
A(Ti) aceleracion del espectro objetivo en el perfodo seleccionado 7); A (7)) 
aceleracion del espectro a escalar en el perfodo 7). Al desarrollar la minimizacion 
del error Ej se llega a la siguiente ecuacion para calcular Fj 

F = El Wj AjTi) ^/AliTj) + A^iTj) 

j nw^AKTd + AjCTd) 

Donde A x ,A y son las ordenadas de los espectros a escalar en sentido X 
(Este-Oeste), Y (Norte-Sur), respectivamente. 
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Figura 7 Modification del espectro escalado por el factor F t y F 1 


Los periodos y factores de peso estan dotados en la tabla 3. 


Tabla 2 

Periodos y Pesos considerados para el escalamiento de sismos. 

Ti ( seg ) 

1.8 

2.0 

3.0 

4.0 

Wi 

0,1 

0,3 

0,3 

0,3 


En caso de que el espectro promedio sea menor que el espectro objetivo 
en el rango de los periodos de interes se multiplica el espectro por un factor de 
escalamiento F ± , en la grafica que se encuentra a la derecha de la figura 7 se 
aprecia el espectro escalado multiplicado por F 1 y se ve que todas las ordenadas 
espectrales son mayores que el espectro objetivo en el rango de interes de los 
periodos. El factor final de escalamiento sera el producto de los dos factores Fj F t . 





Figura 8 Escalamiento de registros slsmicos de estudio al espectro con el sismo 
de diseno DBE. A la izquierda se indica el espectro a escalar y el objetivo; al 
centra el espectro objetivo y el espectro escalado multiplicado por Fj y a la derecha 
el espectro escalado multiplicado por F ; F t . 

3.1 Escalamiento para el sismo DBE 

En la figura 8 se presenta en forma grafica el procedimiento “weighted 
average”; a la izquierda se observa el espectro promedio de los sismos a escalar y 
el espectro objetivo que es el de diseno encontrado con el NEC-15, las lineas 
perpendiculares es el rango de periodos de interes; al centra de la figura 8 se 
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aprecia el espectro escalado multiplicado por Fj se observa que hay ordenadas 
espectrales con valores menores al objetivo por lo que se multiplica por F x . 

A la derecha de la figura 8 se tiene el espectro promedio multiplicado por 
los dos factores y se aprecia que en el rango de interes, todas sus ordenadas 
espectrales son mayores que las del espectro objetivo. En la zona de perfodos 
cortos la diferencia entre estos dos espectros es notable pero esa no es la zona de 
interes, se pudo habertenido valores menores que el espectro objetivo. 

En la tabla 4 se ha utilizado la nomenclatura de EQ1, EQ2, EQ3, EQ4 y 
EQ5 a los sismos de 1966 (Mw 8,1), 1970 (Mw 7,9), 1974 (Mw 8,0), 2007 (Mw 7,9) 
y 1998 (Mw=7.2) En la ultima fila de la tabla se indica el producto de los factores 
Fj F 1 y en la figura 9 se presentan los acelerogramas ya escalados al espectro de 
la Norma Ecuatoriana (NEC-15), con los cuales se encontrara la respuesta en el 
tiempo para el sismo DBE. 


Tabla 4 Calculo de los factores de escalamiento para los sismos de Peru y 
Ecuador con el sismo de diseno DBE. 


FACTORES DE ESCALAMIENTO 


EQ1 

EQ2 

EQ3 

EQ4 

EQ5 

SDNEC15t=1.8 

0.373 

0.373 

0.373 

0.373 

0.373 

SFNormal t— 1.8 

0.027 

0.019 

0.032 

0.315 

0.098 

SFParalel t=1.8 

0.050 

0.022 

0.077 

0.267 

0.161 

SDNEC15t=2 

0.336 

0.336 

0.336 

0.336 

0.336 

SFNormal t-2 

0.032 

0.019 

0.028 

0.370 

0.077 

SFParalel t = 2 

0.041 

0.021 

0.045 

0.242 

0.103 

SDNEC15t=3 

0.224 

0.224 

0.224 

0.224 

0.224 

SFNormal t— 3 

0.023 

0.011 

0.014 

0.173 

0.029 

SFParalel t=3 

0.056 

0.010 

0.023 

0.471 

0.028 

SDNEC15t=4 

0.168 

0.168 

0.168 

0.168 

0.168 

SFNormal t— 4 

0.022 

0.010 

0.014 

0.078 

0.017 

SFParalel t = 4 

0.032 

0.008 

0.023 

0.208 

0.017 

Fj 

4.983 

12.692 

5.880 

0.632 

2.613 

Fi 

1.700 

1.700 

1.600 

1.500 

1.500 

Fi Fj = 

8.471 

21.577 

9.408 

0.948 

3.919 





















163 


Respuesta en el tiempo de Bloque Estructural 3 de UFA-ESPE 



IlM# ll 


Figura 9 Acelerogramas escalados que se utilizaran para encontrar la respuesta 

en el tiempo ante el sismo DBE. 


3.2 Escalamiento para el sismo MCE 

Se precede en forma similar a la indica en el apartado anterior, teniendo 
presente que el espectro del sismo maximo considerado MCE se obtiene 
multiplicando las ordenadas espectrales por 1.4. 

La figura 10 es similar a la figura 8 y en ella se aprecia en forma grafica el 
escalamiento para el nuevo espectro objetivo. En la tabla 5 se resume en calculo 
de los factores de escalamiento. 



Figura 10 Escalamiento de registros sfsmicos en estudio al espectro con el sismo maximo 
considerado MCE. A la izquierda se indica el espectro a escalar y el objetivo; al centra el 
espectro objetivo y el espectro escalado multiplicado por Fj y a la derecha el espectro 
escalado multiplicado porF ; - F 1 . 
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Tabla 5 Calculo de los factores de escalamiento para los sismos de Peru y Ecuador con el 

sismo maximo considerado MCE. 


FACTORES DE ESCALAMIENTO 


EQ1 

EQ2 

EQ3 

EQ4 

EQ5 

SDNEC15t=1.8 

0.373 

0.373 

0.373 

0.373 

0.373 

SFNormal t— 1.8 

0.027 

0.019 

0.032 

0.315 

0.098 

SFParalel t=1.8 

0.050 

0.022 

0.077 

0.267 

0.161 

SDNEC15t=2 

0.336 

0.336 

0.336 

0.336 

0.336 

SFNormal t— 2 

0.032 

0.019 

0.028 

0.370 

0.077 

SFParalel t=2 

0.041 

0.021 

0.045 

0.242 

0.103 

SDNEC15t=3 

0.224 

0.224 

0.224 

0.224 

0.224 

SFNormal t=3 

0.023 

0.011 

0.014 

0.173 

0.029 

SFParalel t = 3 

0.056 

0.010 

0.023 

0.471 

0.028 

SD N EC15t=4 

0.168 

0.168 

0.168 

0.168 

0.168 

SFNormal t— 4 

0.022 

0.010 

0.014 

0.078 

0.017 

SFParalel t=4 

0.032 

0.008 

0.023 

0.208 

0.017 

Fj 

4.983 

12.692 

5.880 

0.632 

2.613 

Fi 

2.000 

2.000 

1.700 

1.700 

1.700 

Fi Fj = 

9.966 

25.385 

9.996 

1.074 

4.442 


4 ESTRUCTURA DE ANALISIS 

En Agosto de 2015 la Universidad de Fuerzas Armadas ESPE, empezo a 
construir en el campus del Valle de los Chillos, los nuevos Bloques Estructurales 
donde funcionaran los Centros de Investigacion y de Post Grado. Una panoramica 
de los edificios que se estan construyendo se presenta en la figura 1 1 , para divisar 
el Bloque Estructural 3 que sera analizado en los siguientes apartados. 



Figura 11 Panoramica del nuevo campo de Investigacion Cientifica y de Post 
Grado de la Universidad de Fuerzas Armadas ESPE. 

En el Bloque Estructural 3 funcionara el Centro de Investigacion de 
Nanotecnologia; es una estructura irregular en elevacion que tiene 4 pisos en un 
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sector (donde esta el numero 3, en la figura 1 1) y 5 pisos en otro sector, como se 
aprecia en la figura 12. 



Figura 12 Vista en elevacion del Bloque Estructural 3, de la UFA-ESPE 




Portico A Portico C Portico E, F 

Figura 13 Porticos en sentido Y de Bloque Estructural 3. 




Portico 12 Portico 13 y 14 

Figura 14 Porticos en sentido X de Bloque Estructural 3. 
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Las columnas del primer piso son de hormigon armado de 80/80 cm, 
sobre cada una de ellas va el aislador de triple pendulo de friccion, como se 
aprecia en las figuras 12 a 14. Sobre el aislador se tiene la losa de aislacion que 
es la losa del primer piso y de este piso al tope las columnas y vigas son de acero. 
En la figura 13 se presentan los porticos en sentido X, y en la figura 14 los porticos 
en sentido Y. 

El aislador con el que se va a construir el Bloque Estructural 3 es el 
FTP8833/1 2- 12/8-6, descrito con detenimiento en Constantinou et al. (2016) y que 
se indica en la figura 1 5. El aislador como tal es circular con un diametro de 91 .44 
cm, pero en su parte superior e inferior tiene una placa cuadrada de un metro de 
longitud, como se aprecia en la figura 16 que corresponde a un aislador que 
estaba siendo ensayado en EPS (Earthquake Protection System) 


36 ' 



Figura 15 Aislador FPT8833/1 2-1 2/8-6, utilizado en la construction del nuevo 
Centro de Investigaciones Cientificas de UFA-ESPE. 

El aislador de la figura 16 esta sin el recubrimiento de goma exterior con el 
que vienen los aisladores de tal manera que se puede visualizar el aislador interior 
que se desliza sobre las placas exteriores. 



Figura 16 Aislador FPT durante el ensayo realizado en EPS. 
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Se habfa indicado que las columnas cuadradas sobre las que se asienta el 
aislador son de hormigon armado de 80 cm de lado y al tener un aislador de 1 .00 
de lado, para que no quede ningun lado en voladizo se construyo una canastilla 
cuadrada de hormigon de 1 .20 m, de longitud que se indica en la figura 1 7. 


1.20 

-+ 4 - 



Figura 17 Canastilla de hormigon de 1 .20 m, de lado sobre columnas cuadradas 
de hormigon armado de 80 cm, de lado. 

Finalmente para acabar con la description de los elementos estructurales 
principales del Bloque Estructural 3, en la figura 18 se presentan las secciones de 
las columnas y vigas. Se destaca que existen vigas secundarias, terciarias y 
cuaternarias de menores dimensiones a la viga “I” indicada en la figura 18 pero 
que no se las considera en el analisis sismico de tal manera que esto contribuye a 
la sobre resistencia de la estructura. 



5 MARCO TEORICO 

La estructura va a trabajar en el rango elastico, solo los aisladores son los 
que van a incursionar en el rango no lineal. Por este motivo en el presente 
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apartado se presenta el modelo de tres fases propuesto por Me Vitty y 
Constantinou (2015) que se utilizo para encontrar la respuesta en el tiempo. 

Para describir el modelo de tres fases, se necesita presentar las siguientes 
ecuaciones: 


R Ueff = R t - hi Para i de 1 a 4. 


(4) 


dl = 


Rj.eff 

Rt 


(5) 


Donde R t es el radio de curvatura; h t es la altura, presentada en la figura 
2; R Ueff es el radio de curvatura efectivo; d\ es la capacidad de desplazamiento. El 
subfndice i, varfa de 1 a 4. A continuacion se indican las 3 fases o regimes del 
modelo. 


• Regimen I 

El desplazamiento se da solamente, en el aislador interior en las placas 
que tienen radios de curvatura R 2 , R 3 , ver figura 1 5 

0 < u < u* 


U* = 2 Oi - M 2 ) R 2 ,eff 


( 6 ) 


F = 


W 

2 R 2,eff 


U + /J. 2 W 


(7) 


Donde u es el desplazamiento lateral del aislador; F es la fuerza lateral 
aplicada; W es el peso que gravita sobre el aislador. A la izquierda de la figura 19 
se aprecia la base rfgida del aislador interior desplazandose en las superficies 2 y 
3; a la derecha se muestra el diagrama de histeresis correspondiente, que 
relaciona la fuerza lateral con el desplazamiento horizontal. 

















Regimen 1 


Desplazamiento 

Figura 19 Desempeno de aislador en Regimen I. 
Fuente: Me Vitty y Constantinou (2015) 
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• Regimen II 

El movimiento de la base rfgida del aislador interior Mega a los topes y 
empieza a desplazarse en las superficies 1 y 4. Normalmente es en este regimen 
en que trabaja un aislador ante un sismo de intensidad moderada y alta, en las 
estructuras. Las ecuaciones se indican a continuacion y en la figura 20 se 
presenta la curva de histeresis correspondiente. 

u* <u < u** 

u** = u* + 2 d{ (8) 

W , 

F = — (u - u*) + ^ W (9) 



Desplazamiento 


Figura 20 Desempefio de aislador en Regimen II. 
Fuente: Me Vitty y Constantinou (2015) 


• Regimen III 

Este regimen se presenta cuando el sismo es extremadamente fuerte y el 
aislador interior se topa con los topes exteriores, en estas condiciones empieza a 
deslizarse la base rfgida del aislador interior en las superficies 2 y 3. Las 
ecuaciones de calculo y el desempeno se indica en la figura 21 . 


u cap = 2 d{ + 2 d* 2 (10) 

F = — — — (u — u”) H — — (u** — u*) + /i x W Ml) 

2 * 2 , 1 // 2R leff 1 

Dependiendo de la fase en que se encuentra el aislador se determina la 
rigidez efectiva y con ella se procede a encontrar la respuesta en el tiempo. No 
solo que se encuentra la rigidez efectiva sino tambien que se halla el factor de 
amortiguamiento equivalente Z eq 
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Desplazamiento 

Figura 21 Desempeno de aislador en Regimen III. 
Fuente: Me Vitty y Constantinou (2015) 


%eq 


Ed 

4 7 tE s 


( 12 ) 


Donde E D es la energia disipada; E s es la energia elastica. Ante el sismo 
DBE, como se vera mas adelante el desplazamiento lateral de la estructura estuvo 
alrededor de los 30 cm, por lo que se verified que el factor de amortiguamiento que 
se obtiene con el modelo de tres fases, sea parecido al que se hallo 
experimentalmente. (EPS, 2015). El factor de amortiguamiento depende de la 
carga vertical aplicada. 


En la figura 22 se presentan valores de amortiguamiento encontrados 
experimentalmente, para una carga vertical sobre el aislador de 250 T., y sometida 
a un desplazamiento lateral de 12” (30.48 cm). Para el Proyecto de la Universidad 
de Fuerzas Armadas ESPE, se compraron 81 aisladores los mismos que fueron 
ensayados con diferentes valores de carga vertical, en concordancia con las que 
van a trabajar. 

En la figura 22 se compara los factores de amortiguamiento encontrados 
experimentalmente con los valores que se hallan con un modelo propuesto por 
Aguiar et al. (2016). 


El factor de amortiguamiento equivalente es el amortiguamiento para el 
primer modo de vibracion pero para encontrar la matriz de amortiguamiento 
considerando todos los modos de vibracion utilizando el modelo de Wilson, se 
necesita definir el factor de amortiguamiento para todos los modos de vibracion 
(Aguiar, 2012). 

Se sabe que en estructuras con aisladores, en los primeros modos trabaja 
el sistema de aislacion y en los modos superiores la estructura que tiene un 
amortiguamiento de 0.05. Con esta hipotesis se obtuvieron los factores de 
amortiguamiento para todos los modos de vibracion, calculando en forma lineal. 

En la figura 22 se tienen valores nominales, maximos y mmimos del factor 
de amortiguamiento, lo propio sucede a nivel de coeficientes de roce. En este 
articulo se esta trabajando con valores nominales o promedio dejando a futuro el 
analisis ante valores mmimos y maximos del coeficiente de roce. 
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Figura 22 Valores de amortiguamiento encontrados experimentalmente y hallados 
con un nuevo modelo propuesto. Fuente: Aguiar et al. (2016) 

6 RESULTADOS 

Para encontrar la respuesta en el tiempo se utilizo el Metodo /? de 
Newmark para multiples grados de libertad. Como se va a encontrar la respuesta 
en cada uno de los porticos del Bloque Estructural 3, para 10 registros sfsmicos, 
se elaboro un programa en MATLAB que permite calcular con una sola corrida 
todos los sismos y exportar los resultados a una tabla de Excel. Por considerarlo 
de interes este programa se lo presenta en el Anexo A. 

Primero se va a presentar la respuesta de desplazamientos, tanto para el 
sismo DBE como para el sismo MCE, en los aisladores y desde el segundo al 
quinto piso. Se recuerda que los aisladores se hallan sobre las columnas del 
primer piso. En el subsuelo los desplazamientos laterales son muy bajos. 

Posteriormente se indica los espectros de pisos, de igual manera para el 
sismo DBE y para el sismo MCE. 

6.1 Desplazamientos encontrados para sismo DBE 

En las figuras 23 y 24 se presentan las respuestas de desplazamiento 
para los porticos en sentido X: A y C, que tienen 4 pisos. En la parte superior 
izquierda se tiene la respuesta del aislador y en las restantes figuras de cada uno 
de los pisos de la superestructura. 
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Portico A 



Figura 23 a) Desplazamientos del aislador, b) desplazamientos del segundo piso, 
c) desplazamiento del tercer piso, d) desplazamientos del cuarto piso. 


Portico C 



Figura 24 a) Desplazamientos del aislador, b) desplazamientos del segundo piso, 
c) desplazamiento del tercer piso, d) desplazamientos del cuarto piso. 

Se destaca que los desplazamientos de los pisos 2 en adelante se 
encuentran en coordenadas absolutas. 

En las figuras 25 y 26 se presentan los desplazamientos de los porticos de 
5 pisos: E, F. Con lo que finaliza las respuestas de los porticos en sentido X. 
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Portico E 




Figura 25 a) Desplazamientos del aislador, b) desplazamientos del segundo piso, 
c) desplazamiento del tercer piso, d) desplazamientos del cuarto piso, e) 
desplazamientos del quinto piso. 

En sentido Y, se tienen tres porticos y son los identificados por los 
numeros: 12, 13 y 14. Las respuestas de desplazamientos de estos porticos, se 
hallan en las figuras 27 a 29. 

En la tabla 6 se presentan los desplazamientos maximos encontrados en 
cada uno de los porticos, en los aisladores y el sismo con el cual se obtuvo ese 
maximo. Los valores maximos de los maximos se hallan con las dos componentes 
del sismo de Pisco de 2007 y son de 35.57 cm (mayor valor positivo) y 46.73 cm 
(mayor valor negativo). 

La capacidad de desplazamiento del aislador FPT8833/1 2-1 2/8-6 es de 54 
cm. De tal manera que van a tener un buen desempeno sfsmico ante estos sismos 
de larga duracion, con dos fases intensas y con varios picos caracterfsticos en la 
zona descendente del espectro. En otras palabras los aisladores van a tener un 
buen desempeno ante estos sismos que tienen un alto potencial destructive. 
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Portico F 



Figura 26 a) Desplazamientos del aislador, b) desplazamientos del segundo piso, 
c) desplazamiento del tercer piso, d) desplazamientos del cuarto piso, e) 
desplazamientos del quinto piso. 

Portico 12 
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Figura 27 a) Desplazamientos del aislador, b) desplazamientos del segundo piso, 
c) desplazamiento del tercer piso, d) desplazamientos del cuarto piso, e) 
desplazamientos del quinto piso. 

Portico 13 




Figura 28 a) Desplazamientos del aislador, b) desplazamientos del segundo piso, 
c) desplazamiento del tercer piso, d) desplazamientos del cuarto piso, e) 
desplazamientos del quinto piso. 


En las tablas 7 a 10 se presentan los desplazamientos laterales maximos 
positivos y negativos, relativos, con respecto al desplazamiento del aislador y se 
aprecia que son valores pequenos. Lo que demuestra que efectivamente la 
superestructura va a trabajar en el rango elastico. Por este motivo se diseno una 
conexion estandar en la union de las columnas de acero tubulares con las vigas 
de acero tipo “I”. 
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Portico 14 



Figura 29 a) Desplazamientos del aislador, b) desplazamientos del segundo piso, 
c) desplazamiento del tercer piso, d) desplazamientos del cuarto piso, e) 
desplazamientos del quinto piso. 


Tabla 6 Valores de desplazamientos maximos: positivos y negativos para sismo 
DBE, en el aislador 


Piso 

Portico 

q.max(+) 

(cm) 

Sismo 

q.max (-) 
(cm) 

Sismo 

Aislacion 

A 

30,93 

7082-07 

38,57 

7036-66 

c 

29,65 

7082-07 

36,34 

7036-66 

E 

30,57 

7082-07 

38,28 

7036-66 

F 

33,96 

7082-07 

43,76 

7081-07 

12 

35,57 

7082-07 

46,73 

7081-07 

13 

33,19 

7082-07 

43,03 

7081-07 

14 

35,09 

7082-07 

43,08 

7081-07 
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Tabla 7 Valores de desplazamientos maximos: positivos y negativos para sismo 


D 

BE en el sec 

undo piso 

Piso 

Portico 

q.max(+) 

(cm) 

Sismo 

q.max (-) 
(cm) 

Sismo 

2 

A 

1,35 

7082-07 

1,62 

7036-66 

C 

1,6 

7082-07 

1,86 

7036-66 

E 

2,24 

7082-07 

2,91 

7036-66 

F 

1,04 

7082-07 

1,33 

7081-07 

12 

1,07 

7082-07 

1,39 

7081-07 

13 

2,36 

7082-07 

3,12 

7081-07 

14 

1,62 

7082-07 

2,07 

7081-07 


Tabla 8 Valores de desplazamientos maximos: positivos y negativos para sismo 
DBE en el tercer piso 


Piso 

Portico 

q.max(+) 

(cm) 

Sismo 

q.max (-) 
(cm) 

Sismo 

3 

A 

2,48 

7082-07 

2,99 

7036-66 

C 

2,92 

7082-07 

3,45 

7051-74 

E 

4,22 

7082-07 

5,44 

7036-66 

F 

1,99 

7082-07 

2,54 

7081-07 

12 

2,01 

7082-07 

2,6 

7081-07 

13 

4,4 

7036-66 

5,8 

7081-07 

14 

3,01 

7082-07 

3,82 

7081-07 


Tabla 9 Valores de desplazamientos maximos: positivos y negativos para sismo 
DBE en el cuarto piso 


Piso 

Portico 

q.max(+) 

(cm) 

Sismo 

q.max (-) 
(cm) 

Sismo 

4 

A 

3,19 

7082-07 

3,87 

7036-66 

c 

3,68 

7082-07 

4,49 

7051-74 

E 

5,52 

7082-07 

7,04 

7036-66 

F 

2,67 

7082-07 

3,39 

7081-07 

12 

2,61 

7082-07 

3,37 

7081-07 

13 

5,76 

7036-66 

7,49 

7081-07 

14 

3,86 

7082-07 

4,89 

7081-07 
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Tabla 10 Valores de desplazamientos maximos: positivos y negativos para sismo 
DBE en el quinto piso 


Piso 

Portico 

q.max(+) 

(cm) 

Sismo 

q.max (-) 
(cm) 

Sismo 

5 

E 

6,14 

7082-07 

7,89 

7051-74 

F 

3,06 

7082-07 

3,87 

7081-07 

12 

2,92 

7082-07 

3,76 

7081-07 

13 

6,51 

7036-66 

8,46 

7081-07 

14 

4,28 

7082-07 

5,41 

7081-07 


En base a los datos de las tablas 6 a 10 en la figura 30 se presenta los 
desplazamientos maximos, en la parte izquierda de la figura 30, tanto positivos (a 
la derecha) como negativos (a la izquierda). A la derecha de la figura 30 se 
presentan las derivas de piso maxi mas encontradas, de igual manera en los dos 
sentidos. 


Desplazamientos Derivas de Piso 








1 1 1 0— 1 1 


-80 

-60 -40 -20 ( 

Desplazamiento 

—•—Maximos Positivos — • — 

3 20 40 

(cm) 

Maximos Negativos 

60 

-0,8 

-0,6 -0,4 -0,2 0 0,2 0,4 

Deriva (%) 

• Maximos Positivos — ©—Maximos Negativos 

0,6 


Figura 30 Desplazamientos y Derivas de Piso maxi mas encontradas en el Bloque 
Estructural 3 de la UFA-ESPE, ante el sismo DBE. 


6.2 Desplazamientos encontrados para sismo MCE. 

Con formato similar a los resultados presentados en el apartado anterior, 
en las figuras 31 a 37 se presentan las respuestas en el tiempo de cada uno de los 
pisos y de cada portico de la estructura del Bloque Estructural 3 de la Universidad 
de Fuerzas Armadas ESPE pero en este caso ante el sismo maximo considerado 
MCE, que es un mega sismo. 

Notese, en la ultima fila de la tabla 5, los factores por los cuales se han 
multiplicado los registros, para realizar el analisis sfsmico, son factores muy altos. 
De tal manera que se va a encontrar la respuesta ante sismos extremadamente 
fuertes y de larga duracion. 
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Portico A 



Figura 31 Desplazamientos de: a) aislador, b) segundo piso, c) tercer piso, d) 
cuarto piso. Para sismo MCE. 


Portico C 



Figura 32 Desplazamientos de: a) aislador, b) segundo piso, c) tercer piso, d) 
cuarto piso, para el sismo MCE 
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Portico E 



Figura 33 Desplazamientos de: a) aislador, b) segundo piso, c) tercer piso, d) 
cuarto piso, e) quinto piso. Para el sismo MCE 



Portico F 
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Figura 34 Desplazamientos de: a) aislador, b) segundo piso, c) tercer piso, d) 
cuarto piso, e) quinto piso. Para el sismo MCE 

Portico 12 



Figura 35 Desplazamientos de: a) aislador, b) segundo piso, c) tercer piso, d) 
cuarto piso, e) quinto piso. Para el sismo MCE 
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Portico 13 


a 



Figura 36 Desplazamientos de: a) aislador, b) segundo piso, c) tercer piso, d) 
cuarto piso, e) quinto piso. Para el sismo MCE 



Portico 14 
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Figura 37 Desplazamientos de: a) aislador, b) segundo piso, c) tercer piso, d) 
cuarto piso, e) quinto piso. Para el sismo MCE 


En la tabla 11 se indican las respuestas maximas de desplazamientos, 
tanto positivos como negativos encontradas en los aisladores ante el sismo MCE. 
Curiosamente todas las respuestas maximas positivas se dan con el sismo de 
2007, y las negativas con el de 1966, registrado en Peru y practicamente los 
aisladores van a trabajar a su maxima capacidad pero resisten, estos mega 
sismos. 

Tabla 11 Valores de desplazamientos maximos: positivos y negativos para sismo 


IV 

CE en los aisladores 

Piso 

Portico 

q.max (+) 
(cm) 

Sismo 

q.max (-) 
(cm) 

Sismo 

Aislacion 

A 

34,69 

7082-07 

51,31 

7036-66 

c 

33,36 

7082-07 

48,73 

7036-66 

E 

35,8 

7082-07 

50,55 

7036-66 

F 

41,57 

7082-07 

53,55 

7036-66 

12 

42,97 

7082-07 

53,92 

7036-66 

13 

40,55 

7082-07 

53,62 

7036-66 

14 - 

41,01 

7082-07 

53,93 

7036-66 


De las tablas 12 a la 15 se indican los desplazamientos maximos relativos, 
con respecto al movimiento horizontal de los aisladores, desde el piso 2 al ultimo. 


Roberto Aguiar, Majorie Arias, Gregory Pachacama y Evelyn Iza 


184 


Tabla 12 Valores de desplazamientos maximos: positivos y negativos para sismo 
MCE en el segundo piso 


Piso 

Portico 

q.max (+) 
(cm) 

Sismo 

q.max (-) 
(cm) 

Sismo 

2 

A 

1,45 

7082-07 

2,07 

7036-66 

c 

1,74 

7082-07 

2,36 

7036-66 

E 

2,59 

7082-07 

3,64 

7036-66 

F 

1,27 

7082-07 

1,62 

7036-66 

12 

1,3 

7082-07 

1,66 

7036-66 

13 

2,94 

7082-07 

3,76 

7036-66 

14a 

1,93 

7082-07 

2,44 

7036-66 


Tabla 13 Valores de desplazamientos maximos: positivos y negativos para sismo 


MCE en tercer piso 


Piso 

Portico 

q.max (+) 
(cm) 

Sismo 

q.max (-) 
(cm) 

Sismo 

3 

A 

2,8 

7082-07 

3,82 

7036-66 

c 

3,35 

7082-07 

4,3 

7036-66 

E 

4,9 

7082-07 

6,82 

7036-66 

F 

2,44 

7082-07 

3,12 

7036-66 

12 

2,45 

7082-07 

3,09 

7036-66 

13 

5,54 

7036-66 

6,97 

7036-66 

14a 

3,62 

7082-07 

4,48 

7036-66 


Tabla 14 Valores de desplazamientos maximos: positivos y negativos para sismo 


MCE en cuarto piso 


Piso 

Portico 

q.max (+) 
(cm) 

Sismo 

q.max (-) 
(cm) 

Sismo 

4 

A 

3,73 

7082-07 

4,93 

7036-66 

c 

4,39 

7082-07 

5,46 

7036-66 

E 

6,49 

7082-07 

8,85 

7036-66 

F 

3,29 

7082-07 

4,18 

7036-66 

12 

3,2 

7082-07 

4,01 

7036-66 

13 

7,27 

7036-66 

8,95 

7036-66 

14a 

4,7 

7036-66 

5,7 

7036-66 
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Tabla 15 Valores de desplazamientos maximos: positivos y negativos para sismo 
MCE en quinto piso 


Piso 

Portico 

q.max (+) 
(cm) 

Sismo 

q.max (-) 
(cm) 

Sismo 

5 

E 

7,27 

7082-07 

9,79 

7036-66 

F 

3,77 

7082-07 

4,79 

7036-66 

12 

3,58 

7082-07 

4,49 

7036-66 

13 

8,23 

7036-66 

9,92 

7036-66 

14a 

5,22 

7082-07 

6,29 

7036-66 


Al igual que en la figura 30, ahora en la figura 38 se presentan los 
desplazamientos maximos de los maximos, a la izquierda y a la derecha 
las derivadas de piso. 


Desplazamientos 



—•—Maximos Positivos —•—Maximos Negativos 


Derivas de Piso 
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— 30— l 

rr 

^ 2Q 


^7 
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5— 

III 


1 



-1 -0,8 -0,6 -0,4 -0,2 0 0,2 0,4 0,6 0,8 

Deriva (%) 

• Maximos Positivos • Maximos Negativos 


Figura 38 Desplazamientos y Derivas de Piso maxi mas encontradas en el Bloque 
Estructural 3 de la UFA-ESPE, ante el sismo DBE. 


Las derivas de piso son valores bajos, menores al 0.8%, que demuestra 
que tendra un buen comportamiento sismico la estructura y los desplazamientos 
laterales maximos superan ligeramente los 60 cm, en el tope del edificio. 


6.3 Espectros de piso para el sismo DBE 

En el analisis no lineal, se hallo las respuestas de desplazamiento, 
velocidad y aceleracion en cada uno de los pisos. La respuesta en desplazamiento 
ha sido presentada en los apartados anteriores. Ahora con la respuesta de 
aceleracion se encontro su espectro, los mismos que se indican en este apartado 
para el sismo DBE y en el siguiente para el sismo MCE. 

De la figura 39 a la 45 se presenta los espectros de piso, encontrados para 
el sismo DBE en cada uno de los pisos y con cada uno de los 10 registros. Se 
empieza colocando los resultados a partir del segundo piso hacia arriba. 
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Portico 12 



Figura 39 a) Aceleracion Espectral piso 2, b) Aceleracion Espectral piso 3, c) 
Aceleracion Espectral piso 4, d) Aceleracion Espectral piso 5. 

Portico 13 



Figura 40 a) Aceleracion Espectral piso 2, b) Aceleracion Espectral piso 3, c) 
Aceleracion Espectral piso 4, d) Aceleracion Espectral piso 5. 


Portico 14 
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Figura 41 a) Aceleracion Espectral piso 2, b) Aceleracion Espectral piso 3, c) 
Aceleracion Espectral piso 4, d) Aceleracion Espectral piso 5. 

Portico A 



Figura 42 a) Aceleracion Espectral piso 2, b) Aceleracion Espectral piso 3, c) 
Aceleracion Espectral piso 4. 



Figura 43 a) Aceleracion Espectral piso 2, b) Aceleracion Espectral piso 3, c) 
Aceleracion Espectral piso 4. 
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Portico F 



Figura 44 a) Aceleracion Espectral piso 2, b) Aceleracion Espectral piso 3, c) 
Aceleracion Espectral piso 4, d) Aceleracion Espectral piso 5. 


Portico E 



Figura 45 a) Aceleracion Espectral piso 2, b) Aceleracion Espectral piso 3, c) 
Aceleracion Espectral piso 4, d) Aceleracion Espectral piso 5. 

En las figuras 39 a 45 se ha colocado el promedio de los espectros 
hallados con los 10 sismos. Ahora interesa ver los espectros en cada piso pero 
con los valores promedios de cada portico; esto se presenta en las figuras 46 a 49. 
Nuevamente se obtiene el espectro promedio con los valores promedios de cada 
piso y de esa forma se tiene un solo valor. 

En la tabla 16 se resume los valores obtenidos. En efecto, se presenta en 
primer lugar los valores maximos del espectro promedio encontrado en cada 
portico y en cada piso. Al final de cada grupo de datos se presenta el valor 
promedio de las ordenadas espectrales de los porticos. 
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Sequndo Piso 



a 


Portico E 

Portico F 

Portico C 

Portico A 

Portico 14 

Portico 13 

Portico 12 

Promedio 


Figura 46 a) Contiene Espectros promedios del segundo piso de todos los 
porticos del Bloque Estructural 3. 


Tercer Piso 



b 


Portico E 

Portico F 

Portico C 

Portico A 

Portico 14 

Portico 13 

Portico 12 

Promedio 


Figura 47 b) Contiene Espectros promedios del tercer piso de todos los porticos 
del Bloque Estructural 3. 
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Cuarto piso 



-Portico E 
Portico F 
-Portico C 
-Portico A 
-Portico 14 
-Portico 13 
-Portico 12 
- Promedio 


Tiempo (s) 


Figura 48 c) Contiene Espectros Promedios del cuarto piso de todos los porticos 
del Bloque Estructural 3. 


Quinto piso 


9 

8 


d 


Portico E 

Portico F 

Portico 14 

Portico 13 

Portico 12 

Promedio 



Figura 49 d) Contiene Espectros del quinto piso de todos los porticos del Bloque 
Estructural 3. 
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Tabla 16 Valores maximos de los promedios de los espectros de piso y valor 


promed 

io de los promedios. Sismo DBE 

Piso 

Portico 

A.max (g) 

Sismo 

2 

A 

3,51 

7039-70 

c 

3,92 

7039-70 

E 

2,35 

7039-70 

F 

5,73 

7039-70 

12 

2,25 

7039-70 

13 

4,47 

7038-70 

14a 

6,9 

7039-70 

Promedio 

4,16 

Piso 

Portico 

A.max (g) 

Sismo 

3 

A 

5,54 

7039-70 

C 

8,42 

7039-70 

E 

4,71 

7039-70 

F 

7,45 

7051-74 

12 

4,73 

7039-70 

13 

7,47 

7051-74 

14a 

8,81 

7039-70 

Promedio 

6,73 

Piso 

Portico 

A.max (g) 

Sismo 

4 

A 

6,62 

7039-70 

c 

11,67 

7039-70 

E 

7,15 

7039-70 

F 

11,05 

7051-74 

12 

6,91 

7039-70 

13 

10,66 

7051-74 

14a 

8,05 

7051-74 

Promedio 

8,87 

Piso 

Portico 

A.max (g) 

Sismo 

5 

E 

8,83 

7039-70 

F 

13,27 

7036-76 

12 

8,27 

7039-70 

13 

10,66 

7051-74 

14a 

9,67 

7051-74 

Promedio 

10,14 


6.4 Espectros de piso para el sismo MCE 

En las figuras 50 a 56 se presentan los espectros de piso, en cada portico, 
en cada piso, para cada uno de los 10 registros y se halla los espectros 
promedios. Esta grafica sirve para saber los valores maximos de los espectros de 
piso y el sismo con el que se encuentran esos maximos. 
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Portico 12 



Figura 50 a) Aceleracion Espectral piso 2, b) Aceleracion Espectral piso 3, c) 
Aceleracion Espectral piso 4, d) Aceleracion Espectral piso 5. Sismo MCE 

Portico 13 



Figura 51 a) Aceleracion Espectral piso 2, b) Aceleracion Espectral piso 3, c) 
Aceleracion Espectral piso 4, d) Aceleracion Espectral piso 5. Sismo MCE. 


Portico 14 




193 


Respuesta en el tiempo de Bloque Estructural 3 de UFA-ESPE 



Figura 52 a) Aceleracion Espectral piso 2, b) Aceleracion Espectral piso 3, c) 
Aceleracion Espectral piso 4, d) Aceleracion Espectral piso 5. Sismo MCE 

Portico A 



Figura 53 a) Aceleracion Espectral piso 2, b) Aceleracion Espectral piso 3, c) 
Aceleracion Espectral piso 4. Sismo MCE. 

Portico C 



Figura 54 a) Aceleracion Espectral piso 2, b) Aceleracion Espectral piso 3, c) 
Aceleracion Espectral piso 4. Sismo MCE 
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Portico F 



Figura 55 a) Aceleracion Espectral piso 2, b) Aceleracion Espectral piso 3, c) 
Aceleracion Espectral piso 4, d) Aceleracion Espectral piso 5. Sismo MCE 


Portico E 



Figura 56 a) Aceleracion Espectral piso 2, b) Aceleracion Espectral piso 3, c) 
Aceleracion Espectral piso 4, d) Aceleracion Espectral piso 5. Sismo MCE 

Ahora, en las figuras 57 a 60 se presentan los espectros medios hallados en 
cada portico y de estos se obtiene un nuevo espectro medio, denominado en las 
graficas como “promedio”. 

Todos los valores maximos de los espectros promedios se indican en la 
tabla 17, en forma similar a la tabla 16 
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Sequndo piso 



Figura 57 a) Contiene Espectros del segundo piso de todos los porticos del 
Bloque Estructural 3. 


Tercer Piso 



b 


Portico E 

Portico F 

Portico C 

Portico A 

Portico 14 

Portico 13 

Portico 12 

Promedio 


Figura 58 b) Contiene Espectros del tercer piso de todos los porticos del Bloque 
Estructural 3. 
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Cuarto piso 



-Portico E 
Portico F 
-Portico C 
-Portico A 
-Portico 14 
-Portico 13 
-Portico 12 
- Promedio 


Tiempo (s) 


Figura 59 c) Contiene Espectros del cuarto piso de todos los porticos del Bloque 
Estructural 3. 


Quinto piso 



d 


Portico E 

Portico F 

Portico 14 

Portico 13 

Portico 12 

Promedio 


6 


Figura 60 d) Contiene Espectros del quinto piso de todos los porticos del Bloque 
Estructural 3. 
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Tabla 17 Valores maximos de los espectros promedios y valores maximos del 
espectro promedio de los promedio. Sismo MCE 


Piso 

Portico 

A.max (g) 

Sismo 

2 

A 

4,212 

7039-70 

c 

4,704 

7039-70 

E 

2,82 

7039-70 

F 

6,876 

7039-70 

12 

2,7 

7039-70 

13 

5,364 

7038-70 

14a 

8,28 

7039-70 

Promedio 

4,99 

Piso 

Portico 

A. max (g) 

Sismo 

3 

A 

6,648 

7039-70 

C 

10,104 

7039-70 

E 

5,652 

7039-70 

F 

8,94 

7051-74 

12 

5,676 

7039-70 

13 

8,964 

7051-74 

14a 

10,572 

7039-70 

Promedio 

E 8,08 

Piso 

Portico 

A. max (g) 

Sismo 

4 

A 

7,944 

7039-70 

c 

14,004 

7039-70 

E 

8,58 

7039-70 

F 

13,26 

7051-74 

12 

8,292 

7039-70 

13 

12,792 

7051-74 

14a 

9,66 

7051-74 

Promedio 

10,65 

Piso 

Portico 

A. max (g) 

Sismo 

5 

E 

10,596 

7039-70 

F 

15,924 

7036-76 

12 

9,924 

7039-70 

13 

12,792 

7051-74 

14a 

11,604 

7051-74 

Promedio 

1 12,17 
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7. COMENTARIOS Y CONCLUSIONES 


Se encontro la respuesta no lineal de cada uno de los porticos del Bloque 
Estructural 3, ante 5 sismos cuya fuente es interplaca tipo thrust, los mismos que 
primeramente fueron escalados al espectro de diseno y espectro maximo 
considerado con que se realizo el diseno estructural, para que sean compatibles 
en el rango de periodos estipulado por el ATC-07 10, en estructuras con 
aisladores sfsmicos. Para el escalamiento se utilizo la tecnica “weighted average”. 

La no linealidad de los aisladores fue representada con el modelo de tres 
fases de Me Vitty y Constantinou (2015). Se encontro las respuestas de 
desplazamientos en el sistema de aislacion y en cada uno de los pisos, ante los 10 
registros considerados y ante el sismo DBE los aisladores tienen la capacidad de 
desplazarse lo que les demanda estos terremotos trabajando en la segunda fase 
del modelo. Ante el sismo MCE los aisladores van a trabajar en la tercera fase 
pero la capacidad de desplazamiento es mayor que la demana. 

Se seleccionaron 4 sismos registrados en Peru y que dejaron mucho dano 
en las estructuras y son los registrados en 1966, 1970, 1974 y 2007. El quinto 
sismo seleccionado fue el de Bahia de Caraquez registrado en Ecuador en 1998. 
Pero como se indico anteriormente fueron escalados (incrementados) sus 
aceleraciones del suelo por factores bastante altos. 

En resumen el Bloque Estructural 3 va a tener un buen desempeno 
estructural ante sismos interplaca tipo thrust que pueden presentarse en la costa 
ecuatoriana. Las derivas de piso maximas que se hallaron fueron menores a 

0.62% para el sismo DBE y 0.78% para el sismo MCE. 
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ANEXO A 

Se realizo un aporte al sistema de computacion CEINCI-LAB con el 
programa newmarklineal2 el cual reporta una tabla con los valores de las 
aceleraciones y desplazamientos en cada piso para cualquier tipo de portico, es 
decir porticos de varios numeros de piso; igualmente el programa presenta una 
tabla con los valores maximos de aceleraciones, desplazamientos y velocidades 
en cada piso, todo esto para el registro del sismo que se desee analizar. 

Los datos de entrada que requiere el programa son los siguientes: 

• [Tabla,Tabla1]=Newmarklineal2(ux,M,C,K,J,dt,beta) 

ux: Vector que contiene los registros del acelerograma en m/s2 

M: Matriz de masas del sistema (T s2/m) 

C: Matriz de amortiguamiento del sistema (T s/m) 

K: Matriz de rigidez del sistema (T/m) 

J: Vector unitario para caso piano. 

dt: incremento de tiempo con el cual se calcula la respuesta. 

beta: Valor de 1/4 para aceleracion constante, 1/6 para aceleracion lineal y 1/8 
escalonada. 

Los datos de salida del programa se describen a continuacion: 

Tabla: empezando de izquierda a derecha imprime una columna con el resumen de 
tiempo, aceleraciones y desplazamientos dependiendo del numero de pisos que tenga el 
portico para un sismo determinado. 

Tablal: imprime el numero de piso, las aceleraciones maximas, los 
desplazamientos maximos y las velocidades maximas de cada piso, es decir que 
escoge el valor absoluto maximo de cada piso dado en la tabla anterior. 

Ejemplo del uso del programa para obtener los resultados en un portico: 

1 ) Se carga el sismo que se desee analizar 
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disp ( 1 Slsmo de 
disp ( 1 Slsmo de 
disp ( " 51 3mo de 
disp ( 1 Slsmo de 
disp ( 1 Slsmo de 
disp ( 1 5isrr.o de 
disp ( 1 5i3ir.o de 
disp ( ' Sl3mo de 
disp E 1 51sro.o de 
disp ( 1 Slsmo de 


Huacho Feru 1966 component e Norte-Sur = l 1 ) 
Huacho Feru 1966 conponente Este-Ceste = 2 1 } 
Chimbote Feru 1970 componente Norte-Sur = 3 1 } 
Chimbote Feru 1970 componente Este-Ceste = 4') 
Lima Feru 1974 componente Norte-Sur = 5 ' ) 

Lima Feru 1974 componente Este-Ceste = 6') 

Feru 2007 componente Norte-Sur = 7 1 > 

Feru 2007 componente Este-Ceste = 3 1 > 

Bahia 1993 componente Norte-Sur = 9'} 

Bahia 1993 componente Este-Ceste = 10 1 > 


sisnio = 'n/d 1 ; %Nombre con el que se generara la hoja de excel 


jjinput ( 1 Increse ei slsmo que desea analizar ') 


if j^=l %Sismo de Huacho Peru 1966 componente Norte-Sur 
load ns_huachoperu 
sisnio = 'nshuachoperu 1 ; 
np= length (ns_huachoperu) ; 
dT=0.02; % Increments de Tiempo 
for i=l:np 

t [ i ) =i ^ dT ; 

uxl ( i) =7 . 723*9 . 3*'ns_huachoperu ( i ) ; 
iie=iie1; 

end 

end 


2) Se ingresa el valor de beta 

3) Se cargan los archivos de la matriz de rigidez lateral y de masas. 
Obtenidos en otro programa. 

4) Se obtienen las propiedades dinamicas, periodos, frecuencias y modos de 
vibracion. Utilizando el programa orden_eig de la galena de programas de 
CEINCI-LAB. 

5) Hallamos la matriz de amortiguamiento con el programa 
amortiguamiento2. 

6) Finalmente llamamos a la funcion newmarklineal2 para que nos reporte 
las tablas antes mencionadas. 
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beta=l/6; k Aceleracion lineal 
% Matriz de rigidez 
load KL142 
% Matriz de Masas 
load M142 

[T,fi,OM]=orden_eig(KL14,M14) ; 

%zeda del hormig6n=Q . 05, zeda del acero=0.03 

zl=[0.2023; 0.1164; 0.03; 0.03; 0.03; 0.031; ^ingresar los valores de zeda para cada modo 
[Cl=amortiguamiento2 (M14 f fi, OM r zl) 

J= [ 1 ; 1 ; 1 ; 1 ; 1 ; 1 1 ; 

[tabla f tablall=Ne™iark:lineal2 (us r M14 r C r KL14 r J r dT r beta) ; 

tabla Walores de la tabla: me reporta tiempo f aceleracione3 de c/pi30 f desplazamientos de c/pi30 
tablal %Me reporta el nmero de piso f aceleracion maxima y desplazamiento maximo 


Con los datos que se obtienen al correr los registros de los 10 sismos en los 7 
porticos que forman el bloque estructural 3 se procedio a realizar el analisis de los 
resultados, obtener las figuras presentadas anteriormente y los valores de 
desplazamientos maximos para cada portico en cada piso, igualmente se obtuvieron los 
espectros de piso. 

Los programas del sistema de computacion CEINCI-LAB se hallan 
publicados en el portal de la Universidad Laica Eloy Alfaro de Manabf, ULEAM, 
dentro de la carpeta de Investigaciones se encuentra la carpeta Centro de 
Investigaciones en Ingenierfa Sfsmica y dentro de esta se halla otra denominada 

CEINCI-LAB 


